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We study a cosmological model in which the boson dark matter gradually condensates into dark 
energy. Negative pressure associated with the condensate yields the accelerated expansion of the 
Universe and the rapid collapse of the smallest scale fluctuations into many black holes, which 
become the seeds of the first galaxies. The cycle of gradual sedimentation and rapid collapse of 
condensate repeats many times and self-regularizes the ratio of dark energy and dark matter to be 
order one. 
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Introduction 

Recent observations including WMAP for cosmic mi- 
crowave background p] and supernovae Type la Hubble 
diagrams[2j independently indicate that some unknown 
dark energy (DE) component dominates the global cos- 
mic energy density and induces the accelerated expansion 
of the Universe. On the other hand, previous many ob- 
servations manifestly indicate that some unknown dark 
matter (DM) dominates the local energy density, whose 
total abundance turns out to be the same order as that 
of DE. 

Amazingly without explicitly specifying dominant con- 
tents of the Universe, the standard ACDM model can ex- 
plain many of the observational data, especially those re- 
lated with the density fluctuations in the linear stage. In 
the non-linear stage however, the standard ACDM model 
cannot explain the very early formation of astronomical 
objects and re-ionization at around z w 20, which is ob- 
servationally required. Some extra mechanism such as 
biasing is inevitable for the galaxy formation, since the 
gravity itself cannot sufficiently compact matter. 

In this letter, we propose a unified model of DE and 
DM in the context of a cosmic phase transition and fur- 
ther the mechanism for the early formation of non-linear 
objects. The key feature of DE would be the volume- 
independent negative pressure, which guarantees the ac- 
celerated cosmic expansion through the Einstein equa- 
tion. The only known matter which shows negative 
pressure would be some classical coherent objects such 
as scalar fields. As an origin of such coherent object, 
we specifically consider the Bose-Einstein condensation 
(BEC) of a boson field with attractive interaction, and 
we can identify the condensation as DE and the exited 
bose gas as DM. The pressure of BEC for ideal gas is 
volume independent and therefore an introduction of at- 
tractive interaction would easily yields negative pressure. 
If we adopt this scenario, the DE dominance means that 
our universe is in almost the ground state described by a 
macroscopic wave function. Thus the largest scale of the 
Universe, as well as the smallest, would turn out to be 
described by quantum mechanics. 

Basic model we construct would have the following 
general features: (a)Bose gas is introduced as DM which 
initially dominates the energy density and the condensate 
of the boson is identified as DE. (b)The condensate has 



negative pressure due to its potential which reflects the 
attractive interaction. For the spatially uniform compo- 
nent of BEC, this negative pressure works as a cosmolog- 
ical constant and guarantees the accelerated expansion. 
(c)The sedimentation of the condensate slowly proceeds 
in the cosmic evolution. (d)When the energy density of 
BEC exceeds some critical value, the condensates rapidly 
collapses into compact boson stars and black holes, which 
work as the standard cold dark matter and also become 
the seeds of galaxies. Simultaneously a new sedimen- 
tation process begins. This cycle of sedimcntation-and- 
collapse repeats many times. These rapid collapses take 
place well after the photon-decoupling stage, and there- 
fore the large scale structure predicted by the ACDM 
model and actually observed pattern in CMB fluctua- 
tions would not strongly be violated. 
cosmic BEC mechanism 

The BEC occurs when the thermal de Broglie length 

1/2 

XdB = (^nh 2 / (mfcT)) exceeds the mean separation 
length of particles n -1 / 3 , where n — N/V is the mean 
number density of the boson and m is the boson mass. 
The critical temperature is given by 
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At the transition point, the chemical potential \l associ- 
ated with the conservation of particle number N within 
the volume V shows singular behavior fj, — > —0. The con- 
densation is possible only for non-relativistic stage T < m 
and the particle number N is conserved. On this stage, 
the cosmic energy density behaves as 
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where no, To are the number density and the tempera- 
ture at some moment in the non-relativistic stage. In 
deriving eq.ffl. we have assumed the adiabatic evolu- 
tion in the sense that the entropy per particle s/n — 

In (e 5 / 2 (mT) 3/2 / (2tt^ 2 ) 3/2 /n) is conserved. The num- 
ber density dependence of the temperature T cx n 2 / 3 is 
the same as Eq. {Q ■ Therefore the condition T < Tq sets 
the upper limit of the boson mass; if we adopt the value 

o 

p now = 9.44 x 10~ 30 g / cm , then m < 2 eV. Moreover, 
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the boson gas dominance condition in the non-relativistic 
stage sets the lower limit of the mass m > 2 eV, provided 
the boson field has the same temperature as that of ra- 
diation. However, if the boson field is not the thermal 
origin, there is no lower limit for the mass. If the above 
conditions hold, the BEC process starts when the boson 
becomes non-relativistic and it continues afterwards. 

Quantum liquid model of cosmic BEC 

In order to describe the dynamics of BEC, the mean- 
field analysis Q using the Gross-Pitaevskii equation is 
usually adopted. This equation has a form of non-linear 
Schrodinger equation, 
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where ip (x, t) is the condensate wave function, V (x) is 
the potential, g = A-nT^a/m, and a is the s-wave scat- 
tering length. If we decompose the wave function as 
ip = \pne %s , and define the velocity as v = hS7S/m, then 
Eq.© reduces to the continuity equation, and the hy- 
drodynamic equation, 



dv 

m 7~ + V 

at 



+ V + gn 



2m^/n 



AV^ 1=0, (4) 



except the last term which is quantum origin (oc Ti ). 
This term can be neglected if the wave number of the 
mode k satisfies k 2 < k 2 = 2mgn/h 2 ; i.e. large scale 
mode. Further, if we choose the attractive interaction 
(g < 0), BEC can be described as fluid with negative 
pressure. This type of argument on BEC yields the de- 
scription by general liquid with the equation of state 
p = -Ap a . 

In this letter, we introduce the following simplest 
model for the BEC: Boson gas is identified as cold DM 
with the equation of state p = 0, and the condensate as 
DE with p = —p. The sedimentation of BEC in the uni- 
form expanding Universe slowly proceeds with the time 
scale r _1 . This setting is very similar to the chaplygin 
gas modelQ with the equation of state p = —A/ p, except 
that in the latter, DE and DM properties are simultane- 
ously included in this single equation of state of a single 
phase. 

Since the Universe is initially extremely uniform, this 
condensation would also be uniform. The energy density 
of the excited bose gas is diluted by the cosmic expan- 
sion, however, that of condensate is not diluted. This 
is because the work supplied to expand the volume V 
to V + dV is —pdV = pdV, which is exactly the nec- 
essary and sufficient amount of energy to produce the 
new condensate in the region dV with the same energy 
density. Therefore eventually the condensate would dom- 
inate the excited gas component and the expansion law 
of the Universe changes from the decelerated expansion 
to the accelerated expansion. 

In the early stage when the bose gas density dominates 
that of condensate, density fluctuations are controlled by 
the dominant component, i.e. the bose gas, and their 



evolution is described by the standard ACDM model. 
However in the later stage when the condensate density 
dominates the bose gas density, the situation drastically 
changes. The linear perturbation equation for the gauge 
invariant quantity, in our case of the equation of state 
p = —p, reduces to 
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where (•••)' = ^7^, and S k = 5p(k)/p. According to 
this, a small scale mode k > 6H 2 rapidly grows, and 
an almost cosmic horizon scale mode k 2 < 6H 2 slowly 
decays, where the comoving wave number k = k/a is 
defined. Moreover, the smaller the fluctuation scale, the 

faster the growing process: 5 k <x exp (tk/Hy Note that 

this rapid collapse is related with the fact that in the gas 
of negative pressure, there is no sound wave (c| < 0). 
The situation now considering is not the (never growing) 
density fluctuations in the de Sitter space in which the 
pressure is a strict constant. 

Let us consider the non-linear stage of the condensate 
collapse. We suppose a uniform spherically distributed 
over-density region of BEC of the radius r and the density 
p. Since the pressure gradient only exists on the surface 
of the sphere, the surface is isotropically compressed to 
form a dense skin. In this process, the skin region with 
the width dr acquires energy 4irr 2 dr (— p) — Aitr 2 drp 
which is exactly the mass of the skin. Therefore the skin 
soon acquires the light velocity. The skin itself has large 
negative pressure in magnitude and therefore self-focuses. 
Since this skin is still located at the same pressure gra- 
dient, it is further compressed and eventually wipe up 
the whole condensate toward the center. This collapsing 
process can be expressed in the evolution equation of the 
skin radius, 
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where = (47r/3) (tq — J p is the time dependent 
mass of the skin, and the right hand side of Eq.JfjJ) is 
the total force acting on the skin A-nr 2 p. All solutions of 
Eq. @ approach to the constant velocity solution 
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Thus the collapse of the condensate is almost the light 
velocity. 

Self gravity of the sphere would further accelerate the 
skin collapse especially in the later stage. The collapse 
would continue until the Heisenberg uncertainty principle 
begins to support the structure, or a black hole is formed, 
or it bounces back outward if the condensate melts at 
the final stage of the collapse. Anyway the collapsed 
condensate forms localized compact objects classified as 
cold dark matter. 

Gradual sedimentation of the condensation — Self Or- 
ganized Criticality 
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We now turn our attention to the global evolution of 
DE/DM in the expanding Universe. The evolution of 
the various energy densities are governed by the set of 
equations, 

P = p c + P a + Pu £r=t = .^* (8) 

P C =Tp g , f)g = -3H P g -Tp g , pl = ~3Hpi, 

where p c , p g , and pi are the energy densities of conden- 
sate, excited boson gas, and the localized energy density 
after the rapid collapse, respectively. 

This set of evolution equations are valid when the con- 
densate does not dominate the energy density: p c < 
p g + pi. Once it dominates p c > p g + pi after the time 
scale r _1 at around z = z c , inhomogeneous components 
of the condensate would rapidly collapse, and some frac- 
tion of p c is transformed into pi. Then the condition 
Pc < Pg + Pi is recovered and the gradual sedimentation 
of the condensate proceeds again following Eq. JSJ during 
the time scale r _1 . This repeated "chase and collapse" 
process by DE(p c ) and DM(p 9 + pi) self-regularizes the 
ratio of them to fix order unity: p c w p g +pi. This kind of 
autonomous dynamics designated as Self Organized Crit- 
icality (SOC) is widely known and observed in nature in 
various places 0- 

In the late stage when the Universe is locked 
in this SOC phase, let us approximate p c — 
p g + pi. Then the Einstein equation d (t) = 
- (4ttG/3) + 3p) a (t) , (a/a) 2 = (8vrG/3) p has the so- 
lution 

a(t)ott 4 / 3 , p(t)ota(ty 3/2 (9) 

which corresponds to the deceleration parameter q — 
-aaa~ 2 = -1/4. 
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FIG. 1: Time evolution of various densities. Green=p c , 
Red=p 9 , and Blue=p;. We set the parameter T = 0.1 in 
the unit of 87rG/3 = 1, and the transition rate from p c to pi 
is set to be 1/3. After 22 cycles of p c dominance, the scale 
factor has evolved about 200 times. DE(p c ) and DM(p 9 + pi) 
self-regularize their ratio to be of order unity. 

Predictions and Observational tests of the model 
We now examine how our model is checked by obser- 
vations. 



(a) Power spectrum of the density fluctuations 
in the linear regime The collapse of the condensate 
proceeds in the smallest scale. This is because the den- 
sity fluctuation is stronger in the smaller scale, and the 
collapse proceeds with almost the light speed. There- 
fore the BEC model does not affect the power spectrum 
in the linear (large scale) regime however it extremely 
enhances the smaller non-linear regime. Especially the 
many black hole formation from DE in the very early 
stage is the most prominent characteristic of our model. 
Such black holes gradually cluster in much larger scale 
along the standard scenario of ACDM model. 

(b) Power spectrum of CMB The condensation 
is supposed to dominate well after the decoupling time 
(z c < 1000 or T -1 > 4x 10 5 year). Therefore the temper- 
ature inhomogeneity generated within the decoupling era 
would not be changed from the standard ACDM model. 
However the integrated Sacks- Wolfe effect, which origi- 
nates from the non-uniform gravitational potential after 
the decoupling time, has a chance to modify the spec- 
trum. In our BEC model, strong non-linear structure 
forms very rapidly in the smallest fluctuation size, which 
is supposed to be well outside of the linear regime. 

(c) Void-wall structure The collapse of the conden- 
sate repeats many times and the smallest scale fluctua- 
tions dominantly collapse in each process. Although in 
the earlier collapsing process the smallest scale fluctua- 
tions (ex. galaxy scale) dominantly forms non-linearity, 
in the later collapsing process fairly larger scale fluc- 
tuations, much larger than the super clusters, have a 
chance to dominantly form non-linearity. Moreover in 
the earlier collapsing process almost spherically symmet- 
ric objects would be preferably formed and therefore the 
non-linearity is almost point-like. On the other hand 
in the later collapsing process, since the scale is large, 
much irregular fluctuations would dominantly collapse 
and therefore the non-linearity can be plain-like. This 
difference comes from the fact that the collapse is trig- 
gered by the pressure which works on surface and not 
by gravity which works on volume. This plain-like non- 
linearity may characterize the largest scale structure of 
the Universe. 

(d) A first star harbors a boson star. In our 

model, once the condensate dominates the comic den- 
sity, vary rapid collapses occur especially in the smallest 
scale of density fluctuations. Associated with this pro- 
cess, the collapsing object can easily fragment into many 
pieces because the pressure is always negative. Therefore 
we can expect the large amount of collapsing objects. For 
bosons, only the Heisenberg uncertainty principle can 
support the structure against the collapse. This is the 
quantum pressure expressed in the last term of Eq.Q). 
This structure is known as the boson star 6]. The size 
R of the object is order of the compton wave length: 
^comton = 2nh/(mc) ~ 2R. This size must be larger 
than the Schwarzschild radius R > 2GM/c 2 for this ob- 
ject not to collapse into a black hole. These equations 
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yield the critical mass for the boson star. 

M cri ticai ~ m 2 vl lm = Mkaup, (10) 

only below which a structure can exist. For example, 
m = 10~ 5 eV yields the critical mass abut the Earth 
mass: M cr iucai = 10 _5 M Q = M®. These compact struc- 
tures may have captured in the process of the first star 
formation. However even this is the case, such seed bo- 
son star and the condensate should be melted into boson 
gas in the high temperature stellar center, leaving no de- 
tectable relic in principle. 

(e) A first galaxy harbors a giant black hole. 
If the fluctuation mass is larger than the critical mass 
M > M cr iti ca i, then the collapse inevitably continues 
until a black hole is formed. In this process, since 
p < 0, dV < 0, no heating is expected thermodynam- 
ically. However the gravitational energy released in this 
collapse would be GM 2 /R, and if this amount of en- 
ergy is used to heat up the condensate, then the tem- 
perature would be, from NT w GM 2 /R and R w GM, 
T w GMm/R w m. Precisely at this point the boson 
becomes relativistic T « m and the particle number no 
longer conserves. Then the chemical potential is not well 
defined and the condensation melts into the ordinary gas 
with positive pressure. Therefore the boson gas stops 
collapse and violently expands outward. In this process, 
some fraction of the condensate would form a central 
black hole and the rest of the condensate would melt 
into the ordinary gas and expand outward. The gravita- 
tional potential of this structure attracts baryon to form 
a cluster around the central black hole. If the size is 
appropriate there forms a galaxy, which harbors a black 
hole in the center and the boson gas and baryon in the 
outskirts. 

What would happen for the expanding boson gas? The 
expansion of once melted boson gas keeps the relation 
T « GMm/R, p w M/R 3 , and therefore T cx p 1 / 3 . 
Since this temperature is still below the critical temper- 
ature, the boson gas would eventually condense and re- 
collapse. Then the condensation melts and the boson 
expands again. Apparently this bounce repeats multi- 
ple times with dissipation until the gas thermalizes com- 
pletely. 

Since the strong non-linearity is formed in the very 
early stage at around z — z c , no extra biasing process 
is necessary in our model. In this sense our model is 
strongly bottom-up type. Early formed stars and galax- 



ies should be the origin of the re-ionization of the Uni- 
verse around z = z c . 

(f) Ejection of matter The above bounces of the 
condensate need not to be isotropic. Especially when the 
total angular momentum is not vanishing, the bounce 
would preferably strong in the direction of the angular 
momentum. This is because the centrifugal force natu- 
rally forms dilute cone region in the direction of angu- 
lar momentum. Therefore we naturally expect, for each 
bounce, that a pair of blobs of the condensate is ejected 
with relativistic velocity parallel to the angular momen- 
tum. The blobs drag baryon which generally emits ra- 
diation. The ejection of blobs is expected to takes place 
repeatedly many times. 

If the ejection period is sufficiently small and the 
aligned blobs are simultaneously visible, then this may 
be observed as two jets emanating from the center of 
the galaxy toward the direction parallel to the angular 
momentum. The essential characteristic of our model is 
that the created jets are all discrete, i.e. sequence of high- 
speed blobs, and not simply a continuous flow of plasma 
gas. Although at present we cannot quantitatively de- 
scribe such highly non-linear processes, we can naturally 
expect the firm positive correlation among the mass of 
the central black hole, the size of the galaxy, and the 
mean separation of the ejected blobs. 

Summary and Further Developments 

We have developed a unified cosmological model of DE 
and DM using the Bose-Einstein condensation of a light 
boson. Fast formation of compact objects and the self 
organized critical Universe are realised. 

In the context of BEC, there are many laboratory ex- 
periments using alkaline atom gas, in which almost ideal 
BEC is realized^}. Especially the collapse dynamics of 
BEC with negative pressure may be useful to de- 

velop our cosmological BEC model. It would be inter- 
esting if the Gross-Pitaevskii equation supplemented by 
the dissipative term, which represents the transition be- 
tween the condensate and excited gas, would lead to the 
quasi-periodic oscillation of BEC collapse. We also have 
to consider the coherence of BEC from various aspects^- 
These issues are now under investigation. 
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